New way to achieve synchronization in spatially extended systems 
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We study the spatio-temporal behavior of simple coupled map lattices with periodic boundary condi- 
tions. The local dynamics is governed by two maps, namely, the sine circle map and the logistic map 
respectively. It is found that even though the spatial behavior is irregular for the regularly coupled 
(nearest neighbor coupling) system, the spatially synchronized (sometime chaotically synchronized) 
as well as periodic solution may be obtained by the introduction of three long range couplings at 
the cost of three nearest neighbor couplings. 
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Recently there have been lot of interests in the small 
world network systems The small world network 

systems lie in between the regular and the random net- 
works where only a small fraction of long range links 
are introduced at the cost of the equal number of reg- 
ular short range links keeping the total number of links 
conserved. The properties of such systems have been 
thoroughly studied very recently because of the fact that 
many biological, technological or social networks fall in 
this category. Intuitively it is expected that the infor- 
mation or source in the small world network will spread 
quickly through out the system because of the presence 
of some long range links. 

On the other hand a number of phase models have been 
proposed over the recent years to describe the dynamical 
behavior of large population of nonlinear oscillators sub- 
ject to a variety of coupling mechanisms |^ S|. A major 
phenomenon that can be observed is the possibility of 
self synchronization among the members of the popula- 
tion. These can represent the fire flies, heart pace maker 
cells, pancreatic beta cells, neurons etc. as well as the 
circuit arrays |l(J . Different cells (units) in different sys- 
tems may be locally governed by some specific rule. For 
example it may be governed by the sine circle map, logis- 
tic map or some differential equations depending on the 
physical or biological system. Furthermore different cells 
may be coupled to each other in various possible ways. 

Even if the local dynamics is regular, the spatial as 
well as the temporal behavior of the regularly coupled 
system may be irregular. It might happen other way also, 
i.e., the dynamics of regularly coupled chaotic oscillators 
may turn out to be regular. We are interested to con- 
fine ourselves in the parameter regime where the spatial 
behavior of a regularly coupled lattice becomes irregu- 
lar. Then the question we ask is whether it is possible to 
control the spatial irregularity by introducing some small 
fraction of long range couplings at the cost of the equal 
number of nearest neighbor couplings keeping the total 
number of couplings conserved. If so, then that would 
be a new way to control spatial irregularity in a spatially 



extended systems. Furthermore, if the solution turns out 
to be synchronized spatially but chaotic temporally, it 
will find good application to the secret communication 
network. We would like to investigate these aspects in 
one dimensional lattice where the local dynamics of in- 
dividual lattice site is governed either by the sine circle 
map or the logistic map |HJ ). The coupling we consider 
here is the simple unidirectional nearest neighbor cou- 
pling. First we discuss the unidirectionally coupled map 
lattice with the local dynamics defined by the sine circle 
map and then by the logistic map. 

Here we consider a one dimensional periodic lattice of 
size N. The local dynamics at any site is governed by 
the following rule: 



x{t + l) = f(x(t)) =x(t)+Q 
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where x is the dynamical variable, t is the time, f2 is 
the frequency and K is the nonlinearity parameter. The 
map given by Eq. (1) is known as the sine circle map 
and very much similar to the dynamical equation in the 
circuits of Josephson Junction arrays. If all the sites in 
the system are independent of each other, the dynam- 
ics of the individual sites will be governed by Eq. (1). 
However, we are interested in the dynamics of the sys- 
tem when all the sites are directly or indirectly coupled 
to each other. The lattice sites may interact with each 
other in a various possible ways. But, here, we consider 
the situation where the dynamical variable at the i th site 
may be governed by the following rule: 

Xi (t + 1) = (1 - e)f(xi(t)) + ef{x l+1 (t)); 1 < i < N, 
and 

x N (t + 1) = (1 - e)f(x N {t)) + ef{ Xl (t)) (2) 

where e is the coupling parameter. The spatio-temporal 
behavior of the system has been studied extensively. |ll]] 
Even though the local dynamics is chaotic, the global 
dynamics of the system may turn out to be chaotic or 
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regular depending on the parameters, initial condition, 
as well as the nature of the coupling. Our intention is 
to start with a random initial condition and fix the pa- 
rameters such that the dynamics of the system becomes 
chaotic both spatially as well as temporally with the reg- 
ular nearest neighbor coupling as described by Eq. (2). 
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FIG. 1. All Xi's at time £=5001 are plotted as a function of 
site index i for a lattice of size 90 where all the sites are cou- 
pled to its nearest neighbor unidirectionally. The parameter 
values are K = \f2, Q, — 0.3 and e = 0.5 respectively. 



For example, we chose the parameter values as e=0.5, 
fl= 0.3 and K = We consider the lattice size, iV=90. 
We allow the system to evolve for 5000 steps before de- 
termining it's state. The spatial behavior of the system 
described by Eq. (2) is shown in Fig. (1) where x^s are 
plotted as a function of i for t=5001. We see that there is 
no regular spatial behavior. The temporal evolution for 
all the sites are also found to be irregular. The temporal 
evolution at the first site is shown in the Fig. (2) as an 
example. 

We introduce the concept of small world network sys- 
tems in the spatially extended system, namely, coupled 
map lattices. Suppose that we have a periodic one dimen- 
sional lattice of N sites. If the lattice sites are coupled 
with its nearest neighbor site unidirectionally, there will 
be N couplings or bonds. We will use the word bond in 
place of coupling for convenience. A fraction of these N 
bonds are replaced by long range couplings. The replace- 
ment of nearest neighbor bond (coupling) by the long 
range bond is made randomly. In the process of replace- 
ment, we make sure that duplicates are not allowed and 
no part of the system becomes isolated. Under such rear- 
rangements, the total number of bonds remain conserved. 
Since the replacement of bonds are done randomly, there 
will be many possible configurations for fixed number of 
bonds being replaced. Therefore, it is necessary to check 
if any of those configurations produces the synchronized 
solutions. 




FIG. 2. This shows the temporal evolution of the dynami- 
cal variable x at site 1 for the same system with same param- 
eters as used in Fig. 1. 



So, we consider the system where p fraction of nearest 
neighbor bonds are replaced by randomly chosen long 
range bonds. The parameter values are are chosen to be 
same as in Fig. (1). Initially all the sites are populated 
randomly and then they evolve according to Eq. (2) with 
the replacement of p fraction of bonds. The solution ap- 
pears to be stable after 5000 iterations. Since there are 
large number of possible configurations, the numerical 
experiment is carried out for 1000 configurations for a 
fixed value of p. The simulation is also carried out for 
different values of p between and 1. It is found that 
only two to ten percentage of configurations lead to syn- 
chronized solutions for each value of p. For p = 0.5 more 
number of configurations produce synchronized solution. 

Thus we experience that there exist some configu- 
rations which may produce the synchronized solutions. 
Therefore, the control of the spatially chaotic behavior is 
possible by the above mentioned mechanism but the per- 
centage value of configurations producing the synchro- 
nized solutions are very low. So, one has to go through 
various sets of re-wiring to achieve the synchronization 
and is therefore difficult to implement in the physical 
systems. One needs to know a re- wiring procedure where 
least number of bonds will be disturbed or replaced. Fur- 
thermore, one should know a systematic way of replace- 
ment of bonds to achieve the synchronization. Intuitively 
we suspect that the replacing of the nearest neighbor 
bonds may be made in a systematic manner and think 
that the bonds should be broken at regular interval in the 
lattice and the disturbed vertices should be connected in 
a regular fashion. More explicitly, let us divide the lattice 
into n number of segments, each segment has length (L) 
equal to — . So, the bond between sites i and i + 1 should 
be broken and the i th site should be connected to (i+L) th 
site. Next, the bond between sites (i + L) and (i + L + 1) 
should be broken and (i + L) th site should be connected 
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to (i + 2L) th site and so on. This procedure should be 
continued until the full lattice is covered. We note that if 
n = N or L = 1, then the modified lattice remains same 
as the regular one. Therefore, we expect that n should 
be between 1 to N — 1. We therefore do the numerical 
experiment for ti=1, 2, 3, 4 and so on. In fact we find 
that n=l and 3 cases work very well. Although, ideally, 
n=l, would be the most suitable solution, we also study 
?i=3, since it is the most ideal case for studying multiple 
rewirings. We find that n=2 and rt=4 have lesser proba- 
bility to synchronize than n=3, over wide regions of the 
SI — e phase space. We will discuss these concepts in de- 
tails in a forthcoming article. For n=3, the re-wiring rule 
is given explicitly for a lattice of size 90: 
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FIG. 3. All Xi's at time t — 5001 are plotted as a function 
of site index i for a system of size 90 where three long range 
couplings are introduced at the cost of three nearest neighbor 
couplings. The parameter values are same as used in Fig.l 



The bonds between sites 1 and 2, 30 and 31, and 60 and 
61 are broken and the links between sites 1 and 30, 30 
and 60, 60 and 90 are established. Furthermore, any ar- 
bitrary site can be chosen as site 1. The spatial behavior 
for this configuration is shown in Fig. (3). The param- 
eters values are same as used for Fig. (1). We notice 
that total synchronization is achieved. The synchroniza- 
tion behavior can be obtained for a lattice of any size 
with the above mentioned rule. The only difference is 
that for larger lattice it takes longer time to achieve the 
synchronization. This is in fact not surprising. However, 
the temporal evolution is chaotic. So, it may be called 
chaotic synchronization. We further note that even if the 
site 1 is coupled with any site around the 30th site with 
small width, synchronization is achieved. Thus a small 
defect in re- wiring does not matter in achieving the syn- 
chronization. 

We clearly note that by disturbing only three regular 
nearest neighbor bonds and establish three long ranged 
bonds in a regular way, one can achieve synchronization 



over a wide range of parameters. The phase diagram in 
the e — SI plane for k = y2 is shown in Fig 4. The shaded 
region (in fig. 4) in the e — SI plane are the allowed pa- 
rameter regime where one achieves the synchronized so- 
lution starting from a random initial conditions. We see 
that there is a large region in the e — SI plane where the 
synchronization is obtained through the re- wiring mecha- 
nism. The synchronized region for only one broken bond 
is larger than that of the 3 bond case by about 20-30 %. 
However, for any other number of rewirings, the allowed 
part of the phase space is much smaller. 
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FIG. 4. Phase diagram in the fl — e plane for the same 
system as described in Fig. 3. Here K — \f2. The shaded 
region corresponds to the allowed parameter space to obtain 
the chaotic synchronization with the re-wiring mechanism. 



Here we consider the same system where the local dy- 
namics is governed by the Logistic map. The logistic map 
has been extensively used to model a number of physical 
as well as biological phenomenon [p"l| . The map is given 
by 

x(t + l) = f(x(t))=nx(t)(l-x(t)) (3) 

where [i is the parameter controling the local dynamics. 
The the dynamics of the variable Xi at the i th site is 
governed by the Eq. (2) with / (x(t)) as given in Eq. 
(3) where t represents the time, and e is the coupling 
strength. In this system also we find that one may ob- 
tain the spatially periodic solution and sometimes syn- 
chronized solution through the re-wiring rule as stated 
earlier. 

In this system, we find that there is a critical value of 
the coupling strength e for each value of /i. For example, 
all Xi S clXG plotted as a function of e for /i=4 in Figs. 
(5) and (6). Fig. 5 is for the system where all the sites 
are regularly coupled to its nearest neighbor. But Fig. 
6 is for the rewired lattice. We clearly see from Fig. 5 
that there is no spatial regularity. On the other hand, 
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for the rewired lattice (Fig. 6) there is a critical value 
of e, say e cr — 0.35 above which the solution is spatially 
periodic of period four. The spatial nature of the solu- 
tion for the rewired lattice is verified in the entire e — fj, 
plane as shown in Fig 7. There are various regions cor- 
responding to various kinds of solutions. The solution is 
spatially periodic with period one or four in the region I, 
the solution is temporally as well as spatially irregular in 
region II, the solution is spatially periodic of period four 
in region III and the solution is chaotically synchronized 
in region IV. 
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FIG. 5. All Xi's at time t — 5001 are plotted as a function of 
e for a system of size 90 where all the lattice sites are coupled 
to its nearest neighbor unidirectionally. Here fi — 4. 
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FIG. 6. All Xi's at time t — 5001 are plotted as a function 
e for a system of size 90 where three long ranged couplings 
are introduced at the cost of three nearest neighbor couplings. 
Here \x = 4. 



Thus we see that the re-wiring mechanism helps us to 
control the spatial irregularity in a large region of the 



parameter space. Furthermore we note the parameter 
region for both the spatial as well as temporal periodic 
solution also large. Therefore, this mechanism with the 
logistic map works well to control spatial and temporal 
irregularity in a large region of parameter space. 
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FIG. 7. Phase diagram in the e — fi plane for the system 
as described in Fig. 6. There are four regions corresponding 
to four kinds of solutions. The solution is periodic of period 
one or four in region I, spatially and temporally irregular in 
region II, spatially periodic of period four in region HI and 
chaotically synchronized in region IV. 



In summary, we have studied spatially extended sys- 
tem where the local dynamics are governed by the sine 
circle map and the logistic map. We have implemented 
the idea of small world network and investigated that 
spatial irregularity in the system can be controlled. We 
have shown that for a regularly coupled system show- 
ing irregular behavior spatially, synchronization can be 
achieved by disturbing a few regular links (nearest neigh- 
bor coupling) and establishing an equal number of long 
ranged links in a regular way. We notice that the re- 
gion in the O — e space where synchronization occurs de- 
pends on how many of these rewirings are done. For one 
and three rewirings, synchronization occurs over large 
regions of the £1 — e space. We will establish the reason 
for this in a forthcoming article. Chaotic synchroniza- 
tion is achieved for the sine circle map with the re- wiring 
mechanism. On the other hand, spatially as well as tem- 
porally periodic solution is achieved for the logistic map. 
Thus we experience that the re-wiring mechanism helps 
us to achieve the chaotic synchronization for some local 
dynamics. The chaotic synchronization is an important 
solution in Josephson junction arrays. Furthermore it is 
very useful as far as the secret communication network 
is concerned. This mechanism may be tested with other 
local dynamics existing in physical or biological systems. 
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